Poisson-noise induced escape from a metastable state 
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We provide a complete solution of the problems of the probability distribution and the escape rate 
in Poisson-noise driven systems. It includes both the exponents and the prefactors. The analysis 
refers to an overdamped particle in a potential well. The results apply for an arbitrary average rate 
of noise pulses, from slow pulse rates, where the noise acts on the system as strongly non-Gaussian, 
to high pulse rates, where the noise acts as effectively Gaussian. 
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I. INTRODUCTION 

Escape from a metastable state underlies a broad range 
of phenomena, chemical reactions, diffusion in solids, and 
population extinction being examples. A classical theory 
of escape was first developed by Kramers [l| . For a Brow- 
nian particle in a potential well, see Fig. [1] he obtained 
the escape rate W in the form W = A esc exp(— AU/ksT) 
and found the prefactor A CS c in a broad parameter range. 
The understanding of the prefactor is important not only 
for completeness of a theory of escape, but also for inter- 
preting the experiment, cf. Refs. 0t3- Therefore much 
effort has been put into extending the Kramers theory, 
see Refs. [H-0] for a review. Most of the obtained results 
refer to generalizations of the Kramers model to other 
types of Brownian motion, for example, to Brownian mo- 
tion in a three-dimensional potential |8| or a periodically 
modulated potential [§-[l5[. 




FIG. 1: A sketch of a metastable potential. The barrier height 
is AU = U(qs) -U{q a ). In the time range t T < t < W" 1 the 
probability current from the metastable state j is independent 
of time and the escape rate W = j [l[ . The probability distri- 
bution for Poisson noise is found by matching the asymptotic 
solutions in regions /, II, III. 

In the last few years escape from a metastable state 
has attracted much interest as a means of detecting non- 
Gaussian noise and potentially determining its statistics 
flfjl - |24j . The proposed noise detectors are continuous sys- 



tems; in the experiment there have been used Josephson 
junctions (2fj| . |23| and nanomechanical resonators pij . It 
is important therefore to have a full theory of escape in- 
duced by non-Gaussian noise, which will include both the 
exponent and the prefactor in the escape rate. 

In this paper we study the probability distribution and 
escape induced by Poisson noise. Such noise is often en- 
countered in photon statistics and in the statistics of cur- 
rent through tunnel junctions. The noise can be far from 
the Gaussian limit, which happens if the noise pulses are 
infrequent, with the rate of the order of the reciprocal re- 
laxation time of the system in the absence of noise t' 1 (a 
more precise condition is specified later) . This parameter 
range was studied in the experiment [24| in particular. In 
the opposite limit of frequent small pulses the noise be- 
comes close to Gaussian. This range is often of interest 
for experiments with Josephson junctions. 

The random motion of the system is very different de- 
pending on whether the system has time to relax between 
the pulses or they are too close in time to change the 
state between them. Respectively, the probability distri- 
butions of the system have very different shapes and the 
expressions for the escape rate have very different struc- 
tures. At the same time, in the both cases the noise has 
the same statistics. One might expect therefore that the 
noise-induced fluctuations can be studied within a single 
approach that would apply for an arbitrary pulse rate. 
Developing such an approach and demonstrating how the 
fluctuations change from the well-understood Gaussian 
limit to the opposite limit of well-separated pulses, along 
with the escape problem, provide the major motivation 
for this paper. 

We model the system by a particle in a potential well, 
with the local minimum of the potential corresponding to 
the metastable state in the absence of noise, see Fig. [T] 
Escape occurs if a sufficiently large outburst of noise 
drives the particle over the potential barrier. We fur- 
ther assume that the particle is overdamped, it has no 
inertia. For infrequent noise pulses the probability distri- 
bution is singular in this case [HI, HBl • As a consequence, 
the Kramers approach cannot be applied and a different 
technique has to be used. 

Our analysis is based on the kinetic equation. We de- 
velop an asymptotic method of solving this equation in 
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the case where the Poisson noise is weak on average, so 
that escape is a rare event, the escape rate W <C t~ 1 . 
This technique applies for an arbitrary relation between 
the two parameters that characterize a Poisson noise, the 
appropriately scaled noise intensity and pulse rate. We 
find the probability distribution near the minimum of the 
potential well (point q a in Fig. [1]) and near the top of the 
potential barrier (point qs) as integral transforms of dif- 
ferent types. The obtained expressions are then matched 
to the distribution in the intermediate range, which is 
found using the WKB-type approximation. This gives 
the full probability distribution, and both the exponent 
and the prefactor in W. 

II. THE MODEL 

The dynamics of an overdamped particle, with coordi- 
nate q, driven by a Poisson noise fp(t) is described by 
the Langevin equation 

q = -U'(q) + fp(t)-(fp(t)}, (1) 
/p(t) 

n 

Here, U(q) is a metastable potential. In the absence of 
noise, the particle has a stable state at the local po- 
tential minimum q a and an unstable stationary state 
at the local maximum qs, see Fig. [TJ The character- 
istic relaxation time is t r = A" 1 , where A a = U"(q a ). 
An important class of systems described by an over- 
damped one-dimensional particle in a metastable poten- 
tial are systems near bifurcation points. Their dynam- 
ics display model-independent features, and their noise- 
induced switching has found applications in various areas 
of physics, see Ref. [2(| for a recent review. In this case 
q is the slow variable, q a is the position of the attrac- 
tor along the g-axis, and qs is the position of the saddle 
point. 

Poisson noise fp{t) is chosen in Eq. ([1]) to be of the 
simplest type, a sequence of unipolar pulses of area \g\ 
that occur at independent instants t n . We assume that 
the sign of g is such that the noise pushes the system 
from q a to qs, i.e., (qs — q a ) 1 9 > 0. The parameter g 
is small, which means that many pulses are required for 
pushing the particle over the barrier, 

(qs - q a )/g > l. 

We note that (qs — q a )/g is not necessarily the largest 
parameter of the theory, i.e., the scaled area of the noise 
pulses, even though it is small, is not necessarily the 
smallest parameter. 

The noise-driven dynamics very strongly depend on the 
average rate v at which the noise pulses are repeated. If 
vt r is not large, the discreteness of the pulses is of pri- 
mary importance, the state of the system changes be- 
tween the pulses when they come at the average rate. 
On the other hand, in the limit of high mean pulse rate, 



vt r ^> 1, the noise is effectively a white Gaussian noise 
with intensity D = vg 2 /2. 

For a high pulse rate, (fp(t)} = vg does not have to 
be small, it can exceed the characteristic dynamical force 

I qs ~ qa 

\/t r even for small \g\. On physical grounds, we 
incorporate the average bias from the noise —q{fp(t)} 
into the potential U(q) and study fluctuations induced 
by the deviations from the mean, i.e., by a zero- mean 
random force fp(t) — (fp). 

There is a similarity between Poisson-noise driven dy- 
namical systems and reaction systems. The latter are 
characterized by a large but finite number of elementary 
units, for example, molecules in a stirred chemical reac- 
tor or individuals in a population. The dynamics is con- 
trolled by reactions between them. The elementary reac- 
tions are short events, which are uncorrclated with each 
other and are characterized by rates; these events present 
a Poisson process. A system may have a metastable state, 
which is approached as a result of the most probable reac- 
tion sequence. However, an unlikely yet possible reaction 
sequence can drive the system far away from this state, 
leading to switching to a different stable state. For one- 
species reaction systems, the prefactor in the escape rate 
can be studied using the Kramers technique [27], [28j] (the 
prefactor in the probability distribution was missed in 
Ref. [p). 

As in reaction systems, in our system the noise leads 
to a finite increment of the system coordinate. However, 
our system is continuous, and in addition it is subject 
to a regular force —U'(q). This leads to a qualitative 
distinction of the escape problem from that in reaction 
systems and to the aforementioned inapplicability of the 
Kramers technique. 

We will study escape using the kinetic equation for the 
probability distribution of the system p = p(q,t). It has 
a form [29[ 

d t p(q,t) = d q [(U'(q) + vg)p(q,t)] 

+v\p(q-g,t)-p(q,t)] (2) 

(for completeness, we give a derivation in Appendix). 

A. A qualitative picture of large fluctuations 

We assume that the noise is weak and that initially 
the system is prepared well inside the potential well in 
Fig. [TJ Over time ~ t r = A~ : the system will approach 
the vicinity of the potential minimum q a and will then 
fluctuate about q a for a long time t 3> t r . Eventually 
there will happen a sequence of noise pulses that will 
cause the system to go over the barrier. In the time 
range ( r Cf< W~ x the probability distribution inside 
the potential well p(q) is quasistationary. It is maximal 
near q a (but generally not exactly at q a , see below). 

Of primary interest to us is the tail of the distribution, 
which is formed by large fluctuations to states far from 
the equilibrium position, (q — q a )/g ^ 1- Such fluctua- 
tions are rare and typically last for time <~ t r , which is 
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the characteristic time of the system dynamics. Because 
a single noise pulse shifts the system by g, the number n 
of noise pulses required for reaching a point q inside the 
potential well in time tf ~ t r is determined by expression 

n - vtf ~ (q - q a )/g. (3) 

Here, vtf is the average number of pulses in time tf] the 
corresponding bias has been incorporated into the poten- 
tial U(q), cf. Eq. (TT]); the difference n — vtf characterizes 
the deviation of the noise from the average. In a rare 
fluctuation |n — vtf\ 3> 1, in agreement with the condi- 
tion (q — q a ) I g ^> 1; cf. Ref. [3fJ where the case vt r < 1 
was outlined. 

The difference between the limits of small and large 
mean pulse rates is that \n — vtf\ may be large or small 
compared to the average number of pulses vt f . The case 
n» vtf corresponds to a comparatively low mean pulse 
rate. The case vtf ^> \n — vtf\ S> 1 corresponds to 
the large mean pulse rate limit, in which the noise is 
essentially Gaussian. This is clear from the expression 
for the probability to have n pulses in time tf, 



(vtf) n exp(-vtf)/nl. 



Indeed, the values of n in the two limiting cases lie, re- 
spectively, on the non-Gaussian tail and in the Gaussian 
part of the distribution P„ . 

For n given by Eq. ([3]), we have \np(q) ~ lnP„. In the 
limits of small and large pulse rate this gives, to leading 
order in (q - q a )/g, 



(q - Qa) , (q - q a ) 



lnp(g) 



In- 

9 gvtf 
(q - Qa) 2 



Vt r *C 



q - q a 



2vgH 



Vt r ~> 



f 



<h 
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»1 (4) 



with tf ~ t r (see Sec. Ill B). Equations (Q| show that the 
tail of the distribution is qualitatively different depending 
on the mean pulse rate v. 

Implicit in the estimate Eq. was the assumption 
that the noise is stronger than the regular force. With 
an appropriate t / this assumption gives the right answer 
not too far from q a where the potential U (q) is parabolic. 
The parabolicity of U(q) thus determines the range of 
applicability of Eq. (|4|). Further away from q a Eq. (|4]) 
still describes the leading-order term in In p(q) for a low 
pulse rate, but becomes inapplicable for a high pulse rate. 

Still further away from q a , where q is outside the po- 
tential well, for V < t < W~ l the distribution p(q) 
corresponds to a quasistationary probability current j, 
which is equal to the escape rate, W = j [lj . One of our 
goals is to find j. This will be done by solving Eq. @ 
separately in the three regions indicated in Fig. [T] and 
matching the solutions. 



III. VICINITY OF THE STABLE STATE 

We first consider region I in Fig. [TJ where the system 
is close to the attractor, \q — q a \ <C \qs — q a \- Here, the 



potential U (q) can be expanded in q — q a keeping only 
the quadratic term, with U'(q) ps A a (g — g a ). This makes 
it possible to find an explicit solution of Eq. ^ . For an 
overdamped particle in a parabolic potential driven by a 
unipolar Poisson noise with random pulse area, the prob- 
ability distribution was found in Ref. (25j . In contrast, 
we consider the case of a constant pulse area, which is rel- 
evant for many physical sources of noise, including elec- 
tron or photon noise; the method [25[ does not apply to 
this case and the distribution is different. The standard 
cikonal approximation often used for white-noise driven 
systems 
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also does not apply. 
We seek the solution of Eq. ([2]) with U' = X a (q — q a ) 
in the form of a Fourier integral, 



p(q) = (27T)- 



dwexp [-iuj{q - q a )} Pa{u). (5) 



Substituting Eq. into Eq. ©, we obtain a linear dif- 
ferential equation for p a {w), with solution 



Pa(uj) = exp 



/A: 



-iguj' 



1 



(6) 



Equation ([6]) gives p a {w) in terms of the exponential in- 
tegral (32|. A constant factor in p a is chosen so that to 
satisfy the normalization condition dqp(q) = 1. 

The maximum of p(q) is located close to the equilib- 
rium position q a . The shape of p(q) near the maximum 
strongly depends on the parameter vt r = v/X a . For a 
high rate of noise pulses, v/X a 3> 1, one can expand the 
exponent in Eq. (|6]) to the quadratic term in u), with the 
result 

p(q) w (A Q /2 7 r J D) 1 / 2 cxp [-\ a {q - q a ) 2 /2D] , 

D = vg 2 /2, \U'(q)/vg\ « 1. (7) 

Equation (0 is the familiar Gaussian distribution near 
the minimum of a potential well for a particle driven by 
white Gaussian noise of intensity vg 2 /2. The inequality 
in Eq. ([7]) is the condition that \wg\ <C 1 for the values 
of lo that give the main contribution to Eq. ([5]); it is 
necessary for the expansion in uj to work. 

The shape of the distribution near the maximum 
is qualitatively different in the opposite limit of small 
mean pulse rate, ^A^ 1 <C (qs — qa)/g, where we use 
qs ~ qa as a typical distance on which the potential 
U(q) becomes essentially nonparabolic. From Eq. (|6|), 
for (q — q a + gv\~ 1 )/g < the integrand in Eq. ([5|) has 
no singularities for Im guj > and exponentially decays 
for Im guj — » oo. Therefore, by closing the integration 
contour over uj in the appropriate halfplane, we obtain 
p(q) = for (q — q a + gvX~ 1 )/g < 0. We note that 
q a — gvX^ 1 is still in the region where U(q) is parabolic 
for small mean pulse rate; this is the equilibrium position 
of the "bare" potential Uo(q), i.e., the potential without 
the noise-induced bias, Uo(q) = U(q) + gvq. In contrast, 
in the Gaussian noise limit, Eq. ([?]). this point is generally 
beyond the range where U (q) is parabolic. 
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The singular behavior of p(q) for small to moder- 
ate v/X a is easy to understand: unipolar noise pulses 
push the system only in the direction of positive [q — 
Qa + 9 iy \~ 1 }/g. If the system has no inertia, its qua- 
sistationary distribution should indeed be zero for (q — 
q a + gvX~ 1 )/g < [U [25||. Because \p a (ui)\ decays as 
aj|-^/A a for large ui, if v < X a distribution p(q) displays 
a power-law divergence, 

p(q)^[(q-q a + 9 ^X- 1 )/gV X "- 1 

for (q — q a + gvX^ 1 )/g — > +0, as is also the case for ex- 
ponentially distributed heights of noise pulses J25| . This 
is very different from the smooth Gaussian peak at the 
distribution maximum for large v/X a , cf. Eq. (f7|). 



A. Distribution tail in the harmonic region of the 
potential 

Of significant interest for us is the tail of the distribu- 
tion, where p(q) ^ 1. It lies for {q — q a )/g ~> 1 and is still 
in the harmonic region of U(q) provided the typical width 
of the distribution peak is small compared to the typical 
distance qs — q a on which U(q) becomes nonparabolic. 
The distribution tail can be obtained from Eqs. (0 and 
((()]) for arbitrary v/\ a . On the tail, integration over oj 
can be done by the method of steepest descent. The ex- 
ternum of the integrand is reached for ui lying on the 
imaginary axis, u> — —ip with p = p(q) given by equation 

H(q,p) = 0, H = u(e P9 -pg-1) -pU'(q). (8) 

In the range considered in this Section U'(q) = X a (q— q a ), 
but the form in which Eq. ([8]) is written allows using the 
expression for H(q,p) where U'(q) is nonlinear in q — q a , 
see below. The result of integration over u> can be put 
into the form 

P(q) = (X a p/2Trd p H) 1/2 cxp[-s(q)}, 

»(?)= f dq'p{q'). (9) 

Function —s(q) is the Legendre transform of the expo- 
nent of p a {u) for ui = —ip(q). The method of steepest 
descent applies if s(q) 3> 1. We note that s{q) can be 
thought of as an action of an auxiliary Hamiltonian sys- 
tem with coordinate q, momentum p, and Hamiltonian 
H(q,p). Equation @ gives the Hamilton- Jacobi equa- 
tion for the auxiliary system, H(q,d q s) = 0. 



1. Applicability of the steepest descent method 

The method of steepest descent applies, i.e., it is suffi- 
cient to keep in the exponent of p a (w) in Eqs. (JS]) and © 
only terms quadratic in (uj + ip) , if the terms of higher 



order in (uj + ip) are small. A term (uj + ip) n enters the 
expansion of lnp a (w) with coefficient K n /n\, 

K n = [d:\np a (io)l = _ w . (10) 

One can check that \K n j g n ~ 2 Ki\ increases with pg 
monotonically from 2/n for pg — > to 1 for pg — > 00. Us- 
ing that K2 = —X^dpH/p [this relation has been used 
in deriving Eq. ^] and taking into account the condition 
s(q) 3> 1, one finds the applicability condition of Eq. © 

\gU"(q)/d p H\ « (pg)" 1 , pU'(q) » \U"(q)\, (II) 

with p(q) given by equation H(q,p) = 0. Formally, 
Eq. (JTT]) was obtained for U"(q) = X a , U'(q) = X a (q-q a ), 
but the form in which it is written makes it applicable 
also outside the range of parabolicity of U (q) , see below. 

B. The limits of almost Gaussian and strongly 
non-Gaussian noise 

The explicit form of the distribution on the tail can be 
obtained in the cases of comparatively large and small 
noise pulse rate. The limit of Gaussian noise corresponds 
to v/X a > (q - q a )/g > 1. We have from Eq. © 
p w 2X a (q~q a )/vg 2 -C 1/g and d p H w vg 2 p/2. The con- 
dition of being on the distribution tail, Eq. (jlip . is met 
for X a (q — q a ) 2 jvg 1 3> 1. For such q, Eq. © coincides 
with Eq. ([7]). The exponent of p(q) coincides also with 
the estimate Eq. (0} for large v/ X a provided the typical 
duration of a fluctuation to q is set equal to tf = t r /2; 
however, this is essentially an artifact, the assumption of 
the noise being much stronger than the regular force in an 
optimal fluctuation does not apply in the Gaussian-noise 
limit, see below. 

In the opposite limit, v/X a <C (q — q a )/g, and in partic- 
ular for vjX a < 1, on the tail of the distribution, where 
(q — q a )/g 3> I, we have from Eq. ((5|) exp(pg) 3> 1. 
Then from Eq. © we obtain for s(q) an expression in 
which the leading-order term coincides with the right- 
hand side of the first equation in Eq. ((4]) if we set 
tj = et r . We note that the prefactor in the distribu- 
tion p(q), Eq. now explicitly depends on q. To the 
leading order it is ~ [2ng(q — g a )] -1 / 2 . The first con- 
dition in Eq. ([TT]) reduces to U'(q)/g > U"(q) and is 
satisfied for (q — q a )/g 3> 1; clearly, the second condition 
in Eq. ([TO is then also satisfied. 



IV. THE REGION AWAY FROM THE 
STATIONARY STATES 

We now consider region 77 in Fig. [TJ where coordinate 
q is far from the both stationary positions of noise-free 
motion, \q a — q\, \q$ — q\ ^> \g\. In the spirit of the WKB 
approximation, we seek the quasi-stationary probability 
distribution p(q) in this range in the eikonal form, 

p(q) = exp[-5(«)], (12) 
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We substitute this expression into Eq. ([2]) for p(q) and 
expand S(q-g) « S(q) -gP + (l/2)y 2 P', where P = S' . 
The expansion is justified if S and P vary smoothly on 
the distance ~ g, even though p(q) does not. Respec- 
tively, we will assume that Pg is not necessarily small; 
however, as will be shown later, g 2 \P'\ <S 1. 

We seek P in the form P w P(°) +PW, with | J^t 1 ^ | < 
|P(°)|, and respectively, S w S , (°)+S' (1) . To find the lead- 
ing order term, P' ', we can disregard the term oc P' in 
the expansion of S; we can also disregard U"p compared 
to U'PWp. Then we obtain from Eq. © 



p(q), = = / p(g)dg, (13) 



p(0) 



where is given by equation H(q,p) = with H de- 
fined by Eq. (JSj) ; note that here we do not assume that 
U'(q) is linear in q. 

The term P (1) is given by a linear equation which fol- 
lows from Eqs. © and ([12")) if we keep linear terms in 
p(i) ; p(0)/ = p / Us i n g tn e relation 

p' =dp/dq=pU"(q)/d p H, (14) 
which follows from Eqs. ©, one finds after some algebra 



P^(q) = (2d p H)- 1 - dp H~U"(q) . (15) 
From Eqs. ([r2])-([T5]) we obtain 

p(q)^C(p/d p H) 1/2 e~ s ^, C = (A a /27r) 1 / 2 . (16) 

Here, s(q) is given by Eqs. © and ©■ The constant 
C is chosen in such a way as to match the probability 
distribution (|9]) in the region where U{q) is parabolic but 

(q - qa)/g > l. 

1. Applicability of the eikonal approximation 

The condition of applicability of the eikonal approxi- 
mation is more complicated than in the simple and well- 
known cases of white Gaussian noise or reaction (birth- 
death) systems, cf. Refs. HH. This is because the 
system does not have one small parameter: even though 
g/ils— Qa) is small, we may still have z^A" 1 ;§> (qs~Qa) / 'g- 
We now show that Eq. (fTTj) provides sufficient applica- 
bility conditions in the whole range of interest. 

It is immediately seen from Eq. (fT4|) that, where 
Eq. (fTTj) holds, we have g 2 |P(°)'| <C 1, as assumed. 
To check that |P^| is small and thus to justify the 
used expansion, one can rewrite Eq. fjl 5|) as P' 1 ) = 
pU"U/2{d p H) 2 with n = d 2 v H - 2p- 1 d p H. It follows 
from condition H(q,p) = that the ratio U/gd p H in- 
creases monotonically from 2/3 for < pg <C 1 to 1 for 
pg — >• oo. Therefore 

P« - pgU"/d p H. 



From this relation one can see that Eq. ([TTj) indeed leads 
to \P^'\ <C |p|,l/<7; in fact, one should separately con- 
sider the case pg < 1 , but the overall analysis is straight- 
forward. 



A. Limiting cases 

The explicit form of the distribution p(q) is easy to 
find in the limit of high rate of noise pulses, v 3> U'(q) / g, 
where the noise is perceived by the system as Gaussian. 
In this case from H = we find p = 2U' /vg 2 , and then 
from Eq. (fj"6|) 



p(q) = (A a /27rP) 1/2 e -^- u ^' D . (17) 

This is the Boltzmann distribution for a particle in a 
potential U(q); the effective temperature D = vg 2 /2 is 
equal to the noise intensity. In the range where U (q) is 
parabolic, Eq. (fT?)) coincides with Eq. ([7]) obtained in a 
different way. 

In the opposite limit of low to moderate rate of noise 
pulses, v <C U'(q)/g, we have 



p(q) ^ [X a /27rgU'(q)} 1/2 cxp[- S (q)}, 



(18) 



where s(q) is given by Eq. (|1 3|) with p w 
g- 1 \\n(U'/vg) + \n\n(U' /vg)} 0. In the range U' = 
A a (<? — q a ) Eq. (fT8|l goes over into the result obtained in 
Sec. Ill in a different way. 

For low pulse rate the value of s(q) is close to what 
follows from the simple expression Eq. (j4} if one es- 
timates the duration of the fluctuation in Eq. Q as 
tf ~ (<j» — q a )/U'(q), which is reasonable (the value of s 
depends on ^ logarithmically). We remind that Eq. ((4]) 
was obtained assuming that the optimal train of the noise 
pulses that bring the system to a remote state q gives 
a much stronger force than the regular force. As seen 
from Eq. (|17[) . this assumption does not apply to the 
case where the mean pulse rate is high and the noise is 
perceived as Gaussian. For Gaussian noise, the optimal 
force is known [33[, for white noise it is equal to twice 
the regular force. 



V. VICINITY OF THE LOCAL POTENTIAL 
MAXIMUM 

Near the local potential maximum, region JiT in Fig.JlJ 
the potential is parabolic, U(q) « U(qs) — \s(g — qs) 2 /2. 
The width of the region where the parabolic approxi- 
mation applies largely exceeds \g\. In a part of this re- 
gion deep inside the potential well the distribution is de- 
scribed by Eq. (|16p . For a small mean rate of noise pulses, 
u/Xs ^ 1, the distribution has an important qualitative 
feature, the square-root divergence of the prefactor for q 
approaching the stationary state qg, which is seen from 
Eq. (TT8"|) . This divergence is smeared out at a distance 
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~ \g\ from qs, see below, but it makes it impossible to 
use the Kramers method of finding the distribution near 
the potential maximum since the method substan- 
tially relies on the smoothness of the distribution on a 
scale that largely exceeds g. 

Sufficiently far outside the well, on the other hand, the 
situation is similar to that discussed by Kramers. The 
effect of weak noise can be disregarded here. The distri- 
bution describes a coordinate-independent quasistation- 
ary probability current j from the potential well. Close 
to qs, but for (q - qs)/g > 1 



W =j ps X s (q - qs)p(q)- 



(19) 



To simplify notations, we assume here that g > 0; the 
generalization to the case g < is straightforward. Equa- 
tion (Unj) follows from the kinetic equation, Eq. ^ , upon 
integration over q from — oo to a given q > qs, with ac- 
count taken of the smoothness of p(q) on a distance ~ g. 
The disregarded correction is oc vg 1 /Xs(q — qs) 2 -C 1; 
note that we do not assume that vg/Xs(q — qs) is small, 
the terms linear in g drop out from the expression for j. 

In order to allow for the singular behavior of p(q) for q 
approaching qs from inside the well, we seek the solution 
of Eq. © with U' = —Xs(q — qs) in the form of a Laplace 
transform, as done previously for periodically modulated 
systems driven by white noise [Hf, [l5| , 



p(q) 



dke- k{ - q - qs) l g p s {k). 



(20) 



Substituting this ansatz into Eq. ([2]) and solving the re- 
sulting first-order equation for ps(k), we find 



p s (k) = Csexp 



(21) 



Deep inside the well, but still in the region where U(q) 
is parabolic, the integral over k in Eq. (|20|) can be cal- 
culated by the steepest descent method. The extremum 
with respect to k lies for k — pg, where p = p(q) is de- 
termined by Eq. © with U' = — Xs(q — Qs)- The result 
of the integration over k is 



p(q) ps C s 



2jr\sg 2 p 



d p H 



1/2 



exp 



- / P(q)dq 



'IS 



(22) 



We note that p(q) increases with q moving inside the well, 
i.e., with increasing (qs — q)/g- 

An analysis completely analogous to that in 
Sec. Ill A 1 shows that the steepest descent method ap- 
plies for sufficiently large (qs — q)/g so that there hold 
inequalities . For (qs — q)/g ^ 1 and for low mean 
pulse rate v/Xs < 1, we have exp(pg) 3> 1, which means 
that the distribution p(q) changes by a large factor when 
q changes by g. The smoothening of the distribution oc- 
curs in the narrow range \(qs — q)/g\ iS 1; the steepest 
descent method does not work in this range. 



Equations (|T6|) and ([22]) match in a broad range of q, 
provided we set 



C s = 



1 



2n\g 



(\ a /\s) 1/2 e-Q, 



Q= / p(q)dq. (23) 



On the other hand, outside the potential well for (q — 
qs)/g 2> 1, the major contribution to the integral over k 
in Eq. (|2Tj) comes from small k < g/(q — qs) -C 1. For 
such k, ps(k) ps Cs, and from Eq. (j2U)) p(q) ss Csg/(q - 
qs)- Using the value of Cs from Eq. (|23j) . one obtains 
from Eq. (|19j) an explicit expression for the escape rate, 



(24) 



W= — (\ a Xs) 1/2 exp(-Q). 
Zir 

This is the central result of the paper. 




FIG. 2: The ratio of the scaled probability density at the 
saddle point to the escape rate Gs, Eq. (1251) . 

It is instructive to look at the ratio of the probability 
distribution to find the particle at the local maximum 
of the potential qs and the escape probability W. We 
normalize p(q) by the characteristic diffusion length in 
the Gaussian-noise limit (vg 2 /nXs) 1 ^ 2 and W by Ag and 
introduce function Gs, 



Gs^^sa 2 /*) 1 ' 2 P(qs)/W. 



(25) 



For the chosen normalization, Gs approaches 1 in the 
limit of large v/Xs- On the other hand, for small 
v/X s we have from Eqs. flU, (HU, ([23]), (|24j) Gs oc 
(v/Xs) 1 t 2 \ \n(v/Xs)\- The overall dependence of Gs on 
v/Xs is shown in Fig. [3J It is seen that Gs approaches 
the Gaussian noise limit comparatively slowly, and that 
it significantly differs from the Gaussian-noise result for 
*VA 5 < 1. 



VI. DISCUSSION 

We have studied the probability distribution and the 
rate of escape from a potential well of an overdamped 
particle driven by Poisson noise. The results cover a 
broad range of the average rate of noise pulses v, from 
small rates, where v is comparable to the reciprocal sys- 
tem relaxation time t~ l , to high rates, v 3> t~ x . Noise- 
induced fluctuations are very different in these limits. For 
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vtr 1 the state of the system changes between succes- 
sive pulses, if they come at the average rate, and thus the 
discreteness of the pulses is essential. For vt r ^> 1 where, 
on average, many pulses occur within the relaxation time, 
the system may perceive the noise as effectively Gaussian. 

The noise is assumed weak, so that the rate of noise- 
induced escape W is small compared to t^ 1 . Escape 
results from a large rare fluctuation, which is an un- 
likely sequence of individual noise pulses, with the over- 
all duration of the fluctuation ~ t r . The condition 
Wt r <C 1 requires that the noise pulse area g be small 
compared to the distance between the stationary states, 
(is ~ Qa)/g 3> 1) which means that many noise pulses 
are needed for escape. However, the developed theory is 
not just an asymptotic theory for g — > 0. This makes it 
significantly different from the theories of escape due to 
Gaussian noise where the noise intensity is the smallest 
parameter, cf. Refs. l3ll l33l or due to reaction random- 
ness in reaction systems, where the reciprocal number of 
particles is the small parameter, cf. Refs. H3, [H. 

In our case, the Gaussian noise limit corresponds to 
vt r 3> (qs — q a )/g, which is incompatible with the limit 
g — > 0. Formally, Poisson noise is characterized by two 
parameters, v and g, and it is the interrelation between 
these parameters and the system parameters that leads 
to the rich pattern of fluctuations. 

In the Gaussian noise limit, our results coincide with 
the well known results for this case. The distribution is 
of the Boltzmann form with temperature given by the 
noise intensity D = vg 2 /2. The escape rate is described 
by the Kramers theory In the Gaussian case the 
time interval between the noise pulses in the most proba- 
ble fluctuation leading to escape varies as the fluctuation 
progresses, keeping the velocity of the system equal to 
U'(q) not too close to the stationary states. 

For vt r -C (qs—Qa)/g, on the other hand, the probabil- 
ity distribution differs qualitatively from the Boltzmann 
distribution, it is singular near the maximum. Away from 
the maximum, along with a non-Boltzmann exponential 
factor it contains a coordinate-dependent prcfactor. The 
latter scales as an inverse square root of the potential 
gradient U'(q), Eq. (jig)) . Here, the time interval between 
the noise pulses in the most probable fluctuation leading 



to escape is such that the velocity largely exceeds U'(q) 
far from the attractor and the saddle point. 

The explicit expression for the probability distribution 
allows one to see the evolution of the distribution de- 
pending on the parameter i/t r , i.e., with the noise vary- 
ing from strongly non-Gaussian to effectively Gaussian. 
We note that the standard WKB-type approximation 
does not apply near the maximum of the distribution 
for vt r < (qs - q a )/g- 

The rate of Poisson-noise induced escape has a sim- 
ple form, Eq. (|24p. It contains an exponential factor 
cxp(— Q) and the prefactor (XaXg) 1 ^ 2 /2tt. The exponent 
Q is given by Eqs. <j8j> and (|23|) . In contrast to escape 
due to white Gaussian noise, for Poisson noise Q is not 
determined by the height of the potential barrier, but 
depends on the actual shape of the potential well. For 
vt r <C (qs — q a )/g it is particularly sensitive to the dis- 
tance between the maximum and minimum of the poten- 
tial [2l], Also, Q depends not on the noise intensity 
vg 2 , but separately on v and g. 

The prefactor in the escape rate is independent of the 
noise parameters and is determined only by the shape of 
the potential near its local minimum and maximum. Re- 
markably, it has the same form as the celebrated Kramers 
expression for the case of an overdamped system driven 
by white Gaussian noise This is unexpected, given 
that the shape of the distribution is generally qualita- 
tively different from the Boltzmann distribution, and the 
analysis used to obtain the rate is very different from the 
Kramers analysis. 

I am grateful to A. Kamenev, M. Khasin, and B. Meer- 
son for useful discussions. This research was supported 
in part by NSF grant CMMI-0900666. 



Appendix A: Kinetic equation for a Poisson-noise 
driven system 

We write the probability density as p(q,t) = (S[q — 
q(t)]), where the averaging is performed over realizations 
of noise fp(t), and q(t) is given by the Langevin equation, 
Eq. (JTJ. Then 



1 f f f t+At 

p(q, t + At) = — J dk(cxp <ik[q- q(t)] + ik J 1 



dhU'(q(tr))+ikugAt\l k [fp\), 



(Al) 



h[fp] = cxp 



-ik 



t+At 



dtif P (ti 



For small At, we can replace U' (q(t{j) with U'(q(t)). Since q(t) is independent of fp(t\) for t\ > t, the averaging 
of Ijc [fp] can be done separately from the first exponential in Eq. (|A1|) . Using the explicit form of the characteristic 
functional for a Poisson noise [34, we obtain 



</ fc [/p]> = cxp [-v At (I - e-^)] » 1 - vAt (1 - e~ ikg ) 



(A2) 
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Wc note also that 

i-( / dk {ikAt [U'(q(t)) + vg] } exp {ik [q - q(t))}) « Aid, {[tf'(g) + vg\ p(q, *)} . (A3) 

In the limit At — >• 0, Eqs. (| Al[) (|A3p immediately give kinetic equation ©. 

I 
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